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Abstract 



We present a consistent renormalization of the top and bottom quark/squark sector of 
the MSSM with complex parameters (cMSSM). Various renormalization schemes are 
Kj^ . defined, analyzed analytically and tested numerically in the decays ^2 — ^ biH^ /W~^ 

0^ I (i = 1,2). No scheme is found that produces numerically acceptable results over all 

^^ ■ the cMSSM parameter space, where problems occur mostly already for real parame- 

(^ . ters. Two schemes are identified that show the most robust behavior. A numerical 

analysis of the four partial stop decay widths is performed in our "preferred" scheme, 
f^ ' "w-b, Ah DR". The full one- loop corrections to the corresponding partial decay widths 

^P ■ are evaluated including hard QED and QCD radiation. We find mostly modest cor- 

rections at the one-loop level. 
i^ ■ 
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1 Introduction 

One of the main tasks of the LHC is to search for Supersymmetry (SUSY) [1] . The Mini- 
mal Supersymmetric Standard Model (MSSM) predicts two scalar partners for all Standard 
Model (SM) fermions as well as fermionic partners to all SM bosons. Of particular interest 
are the scalar partners of the heavy SM quarks, the scalar top quarks, tj (i = 1, 2) and scalar 
bottom quarks bj {j = 1,2) due to their large Yukawa couplings. A scalar top quark tj has 
many possible decay modes, depending on the mass patterns of the SUSY particles. Among 
those decay modes are the decays to a scalar bottom quark, bj, and a charged Higgs boson, 
i7+, or W boson, W~^, 

i,^b,H+ (z,j = l,2), (1) 

U^bjW+ (^,j = l,2). (2) 

If these channels are kinematically allowed they can even be dominant if (most of) the other 
decay modes are kinematically forbidden. Consequently, these processes can constitute a 
large part of the total stop decay width, and, in case of decays to a Higgs boson, they can 
serve as a source of charged Higgs bosons in cascade decays at the LHC. 

For a precise prediction of the partial decay widths corresponding to Eq. ([T]) and Eq. ([21), 
at least the one- loop level contributions have to be taken into account. This in turn re- 
quires a renormalization of the relevant sectors, especially a simultaneous renormalization of 
the top and bottom quark/squark sector. Due to the SU{2)l invariance of the left-handed 
scalar top and bottom quarks, these two sectors cannot be treated independently. Within 
the framework of the MSSM with complex parameters (cMSSM) we analyze various bot- 
tom quark/squark sector renormalization schemes, while we apply a commonly used on-shell 
renormalization scheme for the top quark/squark sector throughout all the investigations. 
Special attention is payed to "perturbativity" , i.e. the loop corrections should not be en- 
hanced by large counterterm contributions resulting from an inappropriate renormalization 
scheme. This turns out to be a constraint that is very difficult to fulfill over the whole 
cMSSM parameter range, where it is especially difficult to achieve this simultaneously for 
small and large values of tan/3. 

Higher-order corrections to scalar fermion decays have been evaluated in various anal- 
ysis over the last decade. The simultaneous renormalization of the top and the bottom 
quark/squark sector was taken into account only in a relatively small subset. In Refs. [21E1 
stop and sbottom decays, including the ones to charged Higgs and SM gauge bosons, have 
been evaluated at 0{as) within the MSSM with real parameters (rMSSM). The numerical in- 
vestigation was restricted to relatively low tan /3 values. These calculations are implemented 
in the program SDECAY [1]. A similar analysis in Ref. [5] included electroweak one-loop cor- 
rections, where again only relatively low tan /3 values were considered. The decays of Higgs 
bosons to scalar fermions, including the charged Higgs decays, at the full one-loop level 
within the rMSSM was presented in Refs. [HlEl, indicating very large one- loop corrections 
for large tan /3. An effective Lagrangian approach in the rMSSM for these types of decays 
was given in Ref. [8], with a numerical analysis for tan/3 = 5. 

The renormalization of the top and bottom quark/squark sector has been analyzed also 
in the context of other calculations in the past. A comparison of different renormalization 
schemes within the rMSSM was performed in Refs. OdU] , focusing on large tan /3. One of the 



renormalization schemes considered therein had been used before within the calculation of the 
two-loop bottom quark/squark contributions to the neutral Higgs boson masses [11] which 
are important for large tan/3 values. Within the cMSSM a renormalization was presented 
in Ref. [12], however without an analysis of its practicability. In Refs. [9l[T3] the top and 
bottom quark/squark sector was renormalized within the cMSSM, but only the QCD part 
needed for the presented calculation was considered. Thus no complete top and bottom 
quark/squark sector renormalization has been performed within the cMSSM. Recently a 
renormalization of nearly all sectors of the rMSSM appeared [H]. In this analysis, however, 
the main focus has been on gauge parameter independence. 

Complex phases, as assumed here in the cMSSM, can be relevant for collider observables 
and possibly extracted from experimental data. Scalar top quark branching ratios at a 
linear collider are discussed in Ref. [15]. Concerning LHC measurements, triple products 
involving the decay of scalar top or bottom quarks are analyzed in Refs. [T6HT9] . Finally, 
rate asymmetries are examined in Ref. [2U]. Depending on assumptions about the LHC 
performance it might be possible to extract information on the phases of Mi, At and Ai, at 
the LHC. 

In this paper we analyze the renormalization of the full top and bottom quark/squark 
sector in the cMSSM. We show analytically (and numerically) why certain renormalization 
schemes fail for specific parts of the parameter space. Finally, we explore the one-loop 
effects for the decays ( 1T1I21) for important parts of the cMSSM parameter space in the favored 
renormalization scheme. We present numerical results showing the size of the one-loop 
corrections, especially including small and large tan/3. The evaluation of the partial decay 
widths of the scalar top quarks are being implemented into the Fortran code FeynHiggs [211 - 
^^. A numerical analysis of all scalar top quark decay modes, involving a renormalization 
of all relevant sectors will be presented elsewhere 



2 The generic structure of the quark/squark sector 

The decay channels (ITP]) are calculated at the full one-loop level (including hard QED 
and QCD radiation). This requires the renormalization of several sectors of the cMSSM as 
discussed below. The sectors not discussed in detail are renormalized as follows: 



The gauge and Higgs sector renormalization has been performed following Ref. 
The gauge boson masses, Myy and Mz, as well as the mass of the charged Higgs boson, 
Mjj, has been defined on-shell while the sine squared of the weak mixing angle, s^, is 
defined via the gauge boson masses, s^ = 1 — ikf^//Vf|. The Z factors for the W boson 
field are also determined within an on-shell scheme while the Z factors of the charged 
Higgs boson field are given by a linear combination of the DR Z factors of the Higgs 
doublets (see Ref. [21]). An additional finite Z factor is introduced to fulfill on-shell 
conditions for the external charged H^ field, tan /3 is defined as DR parameter. 

The Higgs mixing parameter /i has been renormalized via an on-shell (OS) procedure 
for the neutralino and chargino sector 



For the renormalization of the electromagnetic charge we require that the renormalized 



667" vertex in the Thomson hmit is not changed by higher order corrections with respect 
to the corresponding tree-level vertex 



A detailed description of our renormalization of all sectors will be given in Ref. [25]. 

In the following we focus on the top and bottom quark/squark sector. The bilinear part 
of the Lagrangian with top and bottom squark fields, t and b, 



Cr 



t, 



{ili'n)^i{Q~{blK)^iy^^ 



''i/b mass 

contains the stop and sbottom mass matrices M^ and M^, given by 



(3) 



M, 



Ql 



m,. 



Mlc2^iT^ - Q,sl 



rrigXq 



rUqX* 
Ml^ + ml + Mlc^pQqsl 



(4) 



with 



Ag — Aq /i K , 



K = {cot/3, tan/3} for q = {t,h}. 



(5) 



M? and M'^ are the soft SUSY-breaking mass parameters, mq is the mass of the corre- 
sponding quark. Qq and T^ denote the charge and the isospin of g, and Aq is the trilinear 
soft SUSY-breaking parameter. The mass matrix can be diagonalized with the help of a 
unitary transformation Ug, 



D, 



U,-M,-UI 



m 



91 











m| 



U, 



^921 ^<?22 



The scalar quark masses, m^j and tti^j, will always be mass ordered, i.e. mq^ < m 



"92- 



7T7,? 
91,2 



1{ml+ml)+-^i+It!mic,, 



T 



Ml-Ml + Mlc2p{TI-2Qqsl 



4m2|Xg|2 



(6) 



(7) 



The parameter renormalization can be performed as follows. 



Mg ^ M,- 



syiq 



which means that the parameters in the mass matrix M^ are replaced by the renormalized 
parameters and a counterterm. After the expansion 5Wiq contains the counterterm part. 



5yiq,, 

SMq^, 



5Ml^ + 2mq5mq - Mlc^p Qq 5 si + (T^ - Qqsl){c2p SM^ + M| fc 

{A* — fihi) druq + mq{6A* — jjidn — k Sjj) 

(5M 



^2/3 J , 



912 

2 

9h 



5Mq,, = SM? + 2mq6mq + Mlc2p Qq 6sl + QqSl{C2p 5Ml + Ml 5c2p) 



with K given in Eq. ^. 



(9) 
(10) 
(11) 
(12) 



Another possibility for the parameter renormahzation is to start out with the physical 
parameters which corresponds to the replacement: 

U, M, Ut ^ U, M, Ut + U, 6M, Ut = (< J|) + (^< /^|) , (13) 

where 5m?^ and (5m? are the counterterms of the squark masses squared. 6Yg is the counter- 
terno to the squark mixing parameter Yq (which vanishes at tree level, Yg = 0, and corre- 
sponds to the off-diagonal entries in Dg = U^MgUt, see Eq. ([6])). Using Eq. ( TT3|) one can 
express (5Mg by the counterterms 5m?^, Sm^^ and 6Yg. Especially for ^M^^j one yields 

5Mg,, = Ul^Ug,,i6ml - 6ml) + U^^Ug^^Yg + Ug,,U*gJY; . (14) 

In the following the relation given by Eq. flTUl) and Eq. (1141) will be used to express either 
SYq, 6Ag or 6mg by the other counterterms. 

For the field renormahzation the following procedure is applied, 

This yields for the renormalized self-energies 

tg^^ik') = E^k') + lik' - mDiSZg,, + SZiJ - Sml , (16) 

±Uk') = S,,,(A;2) + ^{e - ml)5Zg,, + ^{k' - mJ-J^Z^^ - 6Yg , (17) 

S,-,,(A;2) = S,-,,(fc2) + ^(k' - ml)6Zl^ + ^{k' - ml)6Zg,, - 5Y* , (18) 

tg,,{k^) = T^Uk') + \{k' - ml){5Zg,, + 5ZIJ - bm\ . (19) 

In order to complete the quark/squark sector renormahzation also for the corresponding 
quark (i.e. its mass, rug^ and the quark field, q) renormahzation constants have to be 
introduced: 

mg — )■ TTig + bmg , (20) 

a;^g^(l + ^5Z^^/«)u;^g (21) 



with bviig being the quark mass counterterm and 6Z^ and 6Z^ being the Z factors of the 
left-handed and the right-handed component of the quark field g, respectively. a;=p = |(1=F75) 
are the left- and right-handed projectors, respectively. Then the renormalized self energy. 
Eg, can be decomposed into left /right-handed and scalar left /right-handed parts. Eg and 
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g , respectively, 

tgik) =^u;_Sj'(A;2)+ ^a;+Sf (fc^) + u;-Sf (A:^) + a;+Sf (A;^) , (22) 



^The unitary matrix Uq can be expressed by a mixing angle 9q and a corresponding phase ipq. Then the 
counterterm SYq can be related to the counterterms of the mixing angle and the phase (see Ref. [E]). 



where the components are given by 



t^/^iP) = Sj/^(A;2) + ^{SZJ;/^ + SZi^^*) , (23) 



Sf (fc^) = ^'{k') - -fiSZ!; + 5Zf ) - 6m, , (24) 

Sf (fc^) = Sf (fc^) - ^(5Zf + 5Zf ) - 5m, . (25) 

Note that t^^ik"^) = t^^ik"^)* holds due to CVT invariance. 

3 Field renormalization of the quark/squark sector 

We first discuss the field renormalization of the top and bottom quark/squark sector and 
turn to the parameter renormalization in the next section |H 

The field renormalization, meaning the determination of the Z factors, is done within an 
on-shell scheme for squarks and quarks. We impose equivalent renormalization conditions 
for the top as well as for the bottom quark/squark sector: 

(a) The diagonal Z factors of the squark fields are determined such that the real part of 
the residua of propagators is set to unity. 



Re 



d^,Ak') 



= (^ = 1,2). (26) 



dk^ 
This condition fixes the real parts of the diagonal Z factors to 



Re(5% = -Re- ^^^ 



.. . (^ = 1,2). (27) 



Re above denotes the real part with respect to contributions from the loop integral, 
but leaves the complex couplings unaffected. 

The imaginary parts of the diagonal Z factors are so far undetermined and are set to 
zero, 

Im5%, = (z = l,2). (28) 

This is possible since they do not contain divergences. 

(b) For the non-diagonal Z factors of the squark fields we impose the condition that for 
on-shell squarks no transition from one squark to the other occurs. 



ReE^,,{ml) = 0, (29) 



ReS,-,,(m4) = 0. (30) 

This yields 



SZ, 



ReEU^D-SY, 



K-O 



The counterterin 5Yq is determined in the corresponding parameter renormahzation 
scheme. This means the non-diagonal Z factors of the squark fields do also depend on 
the choice of the parameter renormahzation scheme. 

(c) The quark fields are also defined via an on-shell condition. We impose 

hm l±^Ret,{k)u{k) = , lim u{k)Ret,{k) f"" "^l = , (32) 

where u{k), u{k) are the spinors of the external fields. This yields 

Re^Zf-/^ = -R^{sf'/^(m2) (33) 



'q --- 1 ^q y-q 

+ m, 



+ mj 



L'l^2\ I srR'(^2\ 



s,"K) + s.K 



Sf (mj) + Ef (m^) 



niq 



{\m5Zi - \m5Zf) = zRe {^^.""{mD - ^'."^{mD) = 2Im JReSf (m^)} , (34) 



with S'(/c^) = g^2 ■ Choosing also IraSZ^ = — Im^Z^, the imaginary parts of the 
Z factors can be expressed as 

ImSZ^/^ = ±-^Re {Sf (mj) - Sf (mj)} = ±^ Im JReSf (mj)} . (35) 

Note that the renormahzation condition Eq. fl32l) can only be fully satisfied if the 
corresponding quark mass is defined as on-shell, too. 

The Z factors of the quark fields are not needed for the calculation of the considered 
decay modes of the scalar top quarks (see, however, Ref. [25]). 

4 Parameter renormalization of the top and bottom 
quark/squark sector 

Within the top and bottom quark/squark sector nine real parameters are defined: The 
real soft SUSY-breaking parameters M? , M| and M^ , the complex trilinear couplings At 
and Ab and the top and bottom Yukawa couplings yt and yb which both can be chosen to 
be real. (// and tan/3 as well as the gauge boson masses and the weak mixing angle are 
determined within other sectors, see the beginning of Sect. [2]). Note that the soft SUSY- 
breaking parameter M? is the same in the top as well as in the bottom squark sector due 
to the SU{2)l invariance of the left-handed fields. As in Refs. [9l[13], instead of choosing 
the five quantities M? , M| , M^ and yt, yb the squark masses m? , m? , m? as well as the 
top and bottom quark masses rrit, rrib were taken as independent parameters. 

If a regularization scheme is applied which does not break the symmetries of the model, 
it is sufficient to use counterterms which respects the underlying symmetries. Such countert- 
erms are generated by multiplicative renormalization of parameters and fields of the MSSM. 



The parameter counterterms can be fixed by as many renormalization conditions as inde- 
pendent parameters exist [28]. Concerning the top and bottom quark/squark sector we have 
to set nine renormahzation conditions to define all indepedent parameters. 

For the renormalization of the top quark/squark sector we follow Refs. [9l[T3] but we also 
include electroweak contributions. 

We impose five renormalization conditions, (A)-(E), to fix the parameters of the top 
quark/squark sector: 

(A) The top-quark mass is determined via an on-shell condition, yielding the one-loop 
counterterm Smf. 

6m, = ^Re {m, [S,^(m?) + Sf K)] + [sf (m?) + Sf^K)] } . (36) 

(B), (C) The two top squark masses are also defined on-shell, yielding the real counterterms 

5m| = ReE,-Jm| ) (^ = 1,2). (37) 

(D), (E) Finally, the non-diagonal entry in the matrix of Eq. flT3|) is fixed as 

SY, = lRe{S,^,(m|) + S,-^,(m|)} , (38) 

which corresponds to two seperate conditions as SY, is complex. 

The counterterm of the trilinear coupling 6At is then given via the relation of Eqs. (ITUl) 
and ([H]) as: 

6 A, = -i [f/,-,,t/4(K - 6ml) + Ui^^W.JY: + UlU^JY, - {A, - /x* cot /3) 5m,] 
mt 

+ ((5/i* cot /3 - /i* cot^ /3 5tan/3 ) . (39) 

The definition of 5tan/3 and 5/i is indicated in Sect. |2j 

For the bottom quark/squark sector we are left with four independent parameters which 
are not defined yet. We choose the following four renormalization conditions, (i)-(iv): 

(i) The 62 rnass is defined on-shell: 

5m?^ = ReE,^^(mp. (40) 

(ii)-(iv) These three renormalization conditions are chosen according to the different renor- 
malization conditions listed in Tab. [T] and to the corresponding subsections 14. If - 
14.61 They yield the counterterms (5m{,, 5Ai, and 5Yij where only three of these five 
real counterterms are independent (counting each of the complex counterterms, 
(5/lfe and ^Y^, as two real counterterms). The two dependent counterterms can be 
expressed as a combination of the other ones. 



scheme 


"^^1,2 


nib 


Ab 


Yb 


Sect. 


name 


analogous to the t/t sector: "OS" 


OS 


OS 




OS 


14.11 


RSI 


"mb, Ab DR" 


OS 


DR 


DR 




14.21 


RS2 


"rrtft, n DR" 


OS 


DR 




DR 


14.31 


RS3 


"mfe DR, n OS" 


OS 


DR 




OS 


14.41 


RS4 


"ADR, Ren OS" 


OS 




DR 


Ren: OS 


14.51 


RS5 


"Afe vertex, Re^ OS" 


OS 




vertex 


Ren: OS 


14.61 


RS6 



Table 1: Summary of the six renormahzation schemes for the h/h sector investigated below. 
Blank entries indicate dependent quantities. Ren denotes that only the real part of n 
is renormalized on-shell, while the imaginary part is a dependent parameter. The right- 
most columns indicates the section that contains the detailed description of the respective 
renormahzation and the abbreviated notation used in our analysis. 



Applying these renormahzation conditions fixes the counterterms generated by multiplicative 
renormahzation which fulfill the symmetry relations [28]. 

While the 62 mass is defined on-shell, the fei mass receives a shift due to the radiative 
corrections: 



mr 



"l.OS 



m\ + ( (5m? — ReSr {m\ 



(41) 



The term in parentheses is the shift from m? to the on-shell mass squared. The value 



of m- is derived from the diagonalization of the sbottom mass matrix, see Eq. 

i2 



and 



5m? is defined as a dependent quantity [T0l[29]. m? is the on-shell hi mass squared. In 



Ref. [29] the size of the shift was analyzed while in Ref. [TO] bottom squarks appeared only 
as "internal" particles, i.e. as particles inside the loop diagrams. Concerning the scalar top 
quark decay, Eqs. ([1]) and ([2]), we are now dealing with scalar bottom quarks as "external" 
particles, which are defined as incoming or outgoing particles. These "external" particles 
should fulfill on-shell properties. At this point there are two options to proceed: 



{01) The first option is to use different mass values, ?nj^ and "i^ , for the "internal" and 
the "external" particles, respectively, which can cause problems for charged particles 
as, for instance, scalar bottom quarks (see below). 

{02) The second option is to impose a further renormahzation condition which ensures that 
the 61 mass is on-shell: 



5mi = ReSr (mr 



(42) 



In this case the input has to be chosen such that the symmetry relations are fulfilled 
at the one-loop level. 



As mentioned above, the option (Ol) leads to a problem. The IR-divergences originating 
from the loop diagrams involve the "inner" (i.e. tree-level) mass mi . These have to cancel 
with the real Bremsstrahlung IR-divergences, which are evaluated with the help of the "ex- 
ternal" (i.e. one-loop on-shell) mass m^ , which is inserted into the tree-level diagram (the 
result can, as usual, be expressed with the help of the Soft Bremsstrahlung (SB) factor (5sb: 
A^tree X (^SB, See Ref. [27]). Due to the two different sets of masses the IR-divergences do not 
cancel. One way out would be the use of tree-level masses in all diagrams contributing to the 
part 2Re{A^trcc-^ioop}5 i-6. in all loop diagrams and in the hard and soft Bremsstrahlung 
diagrams. However, this would lead to inconsistencies in the evaluation of the complete loop 
corrected amplitude squared ex: (|A^trecP + 2 Re {A^ tree A^ loop}) due to the different masses 
entering the phase space evaluation. A consistent phase space integration requires the use 
of the same "external" masses for all outgoing particles in all parts of the calculation. 

To circumvent the problem of the non-cancellation of IR-divergences we choose the op- 
tion {02) and impose the further renormalization condition Eq. (H2l) . This requires to choose 
an input that restores the symmetries. Relating (Mg)ii of Eq. @j and (Uq^DgUq)ii with 
Dg of Eq. (ini) yields an expression for the soft SUSY-breaking parameter M? (depending 
on the squark flavor), 

MUq) = |f/,-..r< + \U,j'ml - Mlc^iTl - Q,sl) - m\ (43) 



with q = {t, b}. Requiring the SU{2)i relation to be valid at the one-loop level induces the 

!e also Refs. [2l[3l[30]): 

K. ib) = Ml it) + 6Ml (t) - 6Ml (6) (44) 



following shift in M? (see also Refs 



with 



SM^iq) = \U,,y6mi^ + \U,J'6mi - U.^.U^JY, - U.^.U^JY* - 2m,6m 



+ Ml C2/3 Qq Ssl - {T; -Q,sl) {C2p 6 Ml + Ml Sc2p) . (45) 

In other words, everywhere in the calculation the masses and mixing matrix elements coming 
from the diagonalization of the bottom squark mass matrix, see Eq. ([6]), are used with M? (b) 
including the above shift as in Eq. (14^ . In this way the problems concerning UV- and IR- 
finiteness are avoided. (An exception is the field renormalization of the VT-boson field: In the 
corresponding selfenergies the SU{2)l relation is needed at tree-level to ensure UV-finiteness. 
In this case, tree-level bottom squark masses are used.) 

The various renormalization schemes, following the general choice {02), are summarized 
in Tab. [1] and outlined in detail in the following subsections. 

Comparing with the literature, several of the renormalization schemes (or variants of 
them) have been used to calculate higher-order corrections to squark or Higgs decays. The 
older calculations of the loop corrections have all been performed in the rMSSM. 

• A renormalization scheme employing an "OS" renormalization for nih and Yi, was used 
in Refs. [3l[5] for the calculation of stop and sbottom decays. (The calculation of 
Ref. [3] is also implemented in Ref. [4].) In order to check our implementation given in 
Sect. 14.11 we calculated the decay 61,2 — > tiH~ (see Sect. 15. II for our set-up) and found 
good agreement with Ref. [3]. 



A renormalization scheme similar to the real version of RS2, i.e. "mh, Ai, DR" has 
been employed in Ref. [7] for the calculation of Higgs decays to scalar fermions. In 
the scalar top and the Higgs sector they apply an on-shell scheme (partially based on 
Refs. [311[32]), which differs in some points from our renormalization scheme. 



An on-shell scheme was also used in Ref. [33] (based on Refs. [3Tl|3l]) to evaluate the 

decay /^/V {V = W^,Z). 

In Ref. [35], as a starting point, an on-shell renormalization scheme was used for the 
calculation of the electroweak corrections to r(t2 — ^ ti(f)), {(f) = h,H,A). To improve 
the calculation, the parameters mb, rrit, At and Aj, have also been used as running 
parameters. 



Other "early" papers considered QCD corrections to various scalar quark decays |36f - 
l38] . They mostly employed an on-shell scheme for the quark/squark masses and the 
squark mixing angle 6q, where the counterterm to the mixing angle is S6q oc 6Yq. 

The renormalization scheme "Af, vertex. Re Yb OS" is the complex version of the renor- 
malization used in Refs. [IDIIII] for the 0{ab<ys) corrections to the neutral Higgs boson 
self-energies and thus to the lightest MSSM Higgs boson mass, Mh. 



In the following subsections we define in detail the various renormalization schemes. As 
explained before and indicated in Tab. [T]the two bottom squark masses are renormalized on- 
shell in all the schemes, as in Eqs. ( HOj) and fH2|) . and taking into account the shift of M? (6) 
in Eq. i^^. Within the subsections only the remaining conditions and renormalization con- 
stants are defined explicitly (where 5/i and 5tan/3 are defined within the chargino/neutralino 
sector and the Higgs sector, respectively, in all the different renormalization schemes and are 
not discussed any further). 

4.1 On-shell (RSI) 

This renormalization scheme is analogous to the OS scheme employed for the top quark/squark 
sector. 

(ii) The bottom-quark mass is defined OS, yielding the one-loop counterterm 6mb: 

Smb = ^Re {nib [S,^(m^) + Sf (m,^)] + [Sf (m^) + Sf (m^)] } . (46) 

(iii), (iv) We choose an OS renormalization condition for the non-diagonal entry in the 
matrix of Eq. ( 1T3l) . analogous to the one applied in the top quark/squark sector, 
setting 



SYb = lRe{j:-,Jml) + ^bjrnl)} ■ (47) 
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The conditions (i)-(iv) fix all independent parameters and their respective counterterms. 
Analogous to the calculation of the counterterm of the trilinear coupling At, relating Eq. ( TTOj) 
and Eq. (fT4l) yields the following condition for SAb, 



SAb = — [U-,^U-;j6ml - 6ml) + U~,,UIJY: + U^U-^JY, - {A, - /i* tan/3) 5m,] 

+ {6fi* tan (3 + fi* 6tanf3 ) (48) 

with 5??7,? and 6m3 given in Eqs. f H2|) and (HOj) . respectively. 

4.2 rub DR and A^ DR (RS2) 

(ii) The bottom-quark mass is defined DR, yielding the one-loop counterterm Sm,,: 

6m, = ^Re {m, [S^K^) + S«K)]^^^ + [Ef (m^) + Sf«(m^)]^^J . (49) 

The Idiv terms are the ones proportional to A = 2/e — 7e + log(47r), when using 
dimensional regularization/reduction in D = A — e dimensions; 7e is the Euler 



(iii), (ivl 



constant. 

The complex parameter A, is renormalized DR, 



SA, = — 



m, 

+ 



Ul^M., (ReS^^^(m[)|div - ReS^^^(m|)Uiv 



\ UiUh. {ReJ:iJml)U^ + ReE~,Jml)U 
+ I Ui,Ul^^ (ReS^^^(m?JU. + ReE^^^(m|)Ui. 



^(A - /x* tan/3) Re{m, [T.\;{ml) + Sf (m^)]^.^ 



\SL(^2 



[^rk'^b 



+ Sf^(m2)]^.^}j +5/i*|divtan/3 + /i*(5tan/3 . (50) 



All independent parameters are defined by the conditions (i)-(iv) and the corresponding 
counterterms are determined. Solving Eqs. f lTOl) and flT^ for dY, yields 



5Y, 



I Oil I I 012 I 



Ui.MJ^< 



^< 



+ m,b i Ui^^ U^ {^^b ~ A* 5tan/3 — tan (3 djj) — Ui^_^ U^ {SAb — jj,*Staii/3 — tan /3 Sjj,*) 



y-.^MJ'^l - /"tan/?) - U-,^^UlJAb - /i* tan/3)j 6m, 
where 5m? and 5m? are given in Eqs. (H2|) and (140|) . respectively. 



(51) 
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4.3 rub DR and Yt DR (RS3) 

(ii) The bottom-quark mass is defined DR, yielding the one-loop counterterm Snib. 

6m, = ^-Re {m, [S^K) + Sf K)]^^^ + [Sf K) + Sf«(m^)]^^J . (52) 

(iii), (iv) The complex counterterm 6Yh is determined via a DR renormalization condition, 
setting 

5n = ^Re{s^,,(m[)|div + S^,,K)|div} . (53) 



As in Sect. 14. ![ the renormalization conditions (ii), (iii) and (iv) fix the counterterms 
SrUb and 6Yf,, respectively. Together with the renormalization conditions for Sm? and (5m? 
(see Eq. ( 142|) and Eq. ( 140|) . respectively), 6Ab is given by the linear combination of these 
counterterms as 

SAb = -^ [Ui^U-:jSml - 8ml) + Ui^UlJV: + UlU^JY, - (^ - / tan/3) <5m,] 

+ ((5/i* tan /3 + /i* (5tan/3 ) , (54) 

which, of course, shows the same analytical dependence of the independent counterterms as 
5Ab in Eq. (gHD in Sect. O 



4.4 mb DR and yb on-shell (RS4) 

(ii) The bottom-quark mass is defined DR, yielding the one-loop counterterm dm,,: 

6m, = ^Re {m, [^^(m^) + Sf (m^)]^^^ + [Sf (m,^) + Sf«(m^)]^^J . (55) 

(iii), (iv) The complex counterterm 6Yb is fixed by an on-shell renormalization condition, 
as in Sect. 14. ![ 

5Y, = ^Re{^,Jml) + J:,Jml)} . (56) 

As in Sect. 14.11 and in Sect. 14.31 the renormalization conditions (i)-(iv) fix the counterterms 
(5m? , Smb and 6Yb. The further renormalization condition Eq. (H2|) determines the countert- 
erm (5m? . Analogous to Sect. 14.11 and to Sect. 14. 3[ SAb can be expressed in terms of these 
counterterms, 

SAb = -^ [Ub,U-;jSml - Sml) + U-.^U-^JYC + U^UiJYb - {A, - /.* tan/3) dm,] 

+ (5/i* tan /3 + /i* (5tan/3 ) , (57) 

which, of course, has the same form as in Eqs. (148|) and (l54l) . 
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4.5 Ab DR and Re Yb on-shell (RS5) 

(ii) In the subsections 14. H - l4!4] the second renormalization condition defines the bottom 
quark mass. In this scheme, we choose an on-shell renormalization condition for 
the real part of the counterterm 6Yh which determines Re SYf, as following 



Re 5Y, = ^ Re { ReE^^^ (m[ ) + ReS^^^ (mfj } . (5^ 



(iii), (iv) The complex Af, parameter is defined DR 

1 



SAh= — 
nib 

1 . / 



^fcn^L (^^^feii("^fei)ldiv - ReS^^^(m^j|div 



b2' 



+ \UI.M,. (ReS^,,K)|di. + R^S^^^K^)|div)* 
- ^(A - /i* tan/3) Re{m, [S^K^) + T.^{ml)\^^ 

+ [Sf^(m2) + Sf^(m^)]^.J + Vldivtan/3 + /i* 5tan/3 . (59) 

With the conditions (i)-(iv) the independent counterterms 5m? , Re(5Yj, and 6Ab are 
determined, and (5?n? is given by Eq. P2|) . The missing counterterms 5mb and Im^yfe can 
be expressed by the independent counterterms. Relating Eq. (fTOj) . here explicitly written as 

(5M^)i2 = (A^ — /itan/3) 5m;, + mf, (5A^ — /i5tan/3 — (5/itan/3) , (60) 

and Eq. (TT4|) . here with ^Yj, explicitly split into a real and an imaginary part 

{5yi,),, = UlU,^^{5ml-5ml) 

+ Ul^^ U-,^^ (Re 5Yb + ^ Im 5^) + U-,^^ U^^^ (Re 5^ - i Im 5^) , (61) 

results in the two equations 

Re { A^ - ^ tan /?} 5mb = -m^ Re 6Ab - Re 55 - Im | U-*^^ U^^^ - Uj,^^ U^^^ } Im 5^ , (62) 
Im { A^ - /i tan /3} 5mb = +mb Im 5^^ - Im 55 + Re | U-*^^ U^^^ - U^^^ U-*^^ \ Im 5^ (63) 

with 

5S = —rrib (/i5tan/3 + 5/i tan /3) — U^ Ui^^{5m^ — Sm^ ) 

-(uiU,^.+U,^^U-;)Re6Yb, (64) 



where 5?Ti? and 6m3 are given by Eq. P2|) and Eq. fHOl) . 

The above two equations, fl62|) and fl63l) . can be solved for Im5y6 and 5m6, yielding 
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dnib = _ — — , (65) 

H/'pL/'i H/'tL/'p 

lra5Y, = ^^^^^ (66) 



with 



k 


= -Re\W Ui -Ui W 

(_ bii 022 012 b21 


Cr 


= +mb Re 6 Ah + Re6S , 


Ci 


= —nib Im 5Ab + \m.5S . 



ttr = Re {Al — fi tan (3} , (67) 

Oj = Im{y4^ — /itan/3} , (68) 

br = +^^{uiU,.,-U,^uQ , (69) 

(70) 

(71) 

(72) 

4.6 At via vertex and Re Yi, on-shell (RS6) 

(ii) An on-shell renormalization condition is imposed for the real part of the coun- 
terterm 6Yb which determines Re 6Yb as 

ReSYb = ^Re{ReS^^^(m[) + R^S^^^(my} . (73) 

(iii), (iv) The renormalization conditions introduced here are analogous to the prescriptions 
used in Refs. [9HTT], but extended to the complex MSSM. The complex parameter 
Ab is renormalized via the vertex Ab\b2, denoting the renormalized vertex as 
Hp%p1^pI)^ see Fig. m 



h 



A [ j = tA{p%pl,pl] 




Figure 1: The renormalized vertex A(p^,p? ,p? 



fei'-^ba' 
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The tree-level vertex Ab\b2, denoted as V^rt,; , is given as 



Ab[b2' 



V.rn 



terrib 



^^1*2 2Mws„cosl3l 



Ui,U-;jfi cos P,, + a; sin p,,) 
~^h, Ul^, (Ai* cos /3n + Ab sin /3n) 



(74) 



where /3n is the mixing angle of the CP-odd Higgs boson fields with /3n = /3 at 
tree-level. Note that in our renormalization prescription we do not renormalize 
the mixing angles but only tan /3 appearing in the Lagrangian before the transfor- 
mation of the CP-odd Higgs boson fields into mass eigenstate fields is performed. 
The renormalized vertex reads, 



Aip%pl,pQ = A{p\,pl,pQ + 



tenib 



2MwS 



W^w 



tan/3 



Ui.u~;jAi-u,^^uiMb 



+ 



Ui..u~;jf^-u,^^u^^jf^* 



+ 



Ui^Ul, (/^ + tan/^A*) - U-b^Ul^ (/i* + tan/3^ 

^^^ + ^{6Zl^ + 6Zi^^ + 6Zaa) + sin/3 cos/3(5tan/3 



rrib 



+ 



6Z, 



5Ml, 5s, 



'w 



2M^ Sw 



U',^U~;jfi + tan/3A*) - Uj,^U^^Jfi* + tan/3A 

+ tlm{Ui^Ul^ifi + tanf3Al)}6Zj,^^+tlm{Uj,^Ul^ifi + tan(3Al)}6Z-^ 

1 
~ 2 



U-^rn {A; - /itan/3) - U-,m {Ab - /inan/3) 



sz 



AG 



621 

(75) 



The off-diagonal Z factors are determined according to Eq. (I3T 



^\. = +2 ■ 



SZ-; 



621 



-2 



ReS^^^ (m| )-Re6Yb-i Im 6Yb 
(m? — m? ) 

^ bi 62 

ReE^^^ (m?^ )-Re6Yb + i Im 6Yb 
(m? — m? ) 

^ bi 62 ' 



--■.SZ^^ - 



2iIm6Yb 



bi2 ^2 _ j^2 



--■.5Z1 - 



b2 

2i Im 6Yb 



b2l ^ _ J72? 
fei 62 



Introducing appropriate abbreviations we get 

HpI.pI^pI) = HpI,p1.pI) + 

ie 



2Mws, 



{m, tan/3 {Ui^U-;jAl-Ui^^U-;jAb) + 6M + zUY lm6Yb 



+ 



Ub,U-;^, (/" + tan/3A*) - t/j^^f/^^^ (/i* + tan/3A) (5m, + 5Zd) [ + 6Z0 



(76) 



(77) 



with 
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5M = rrib 



UwU^,Jf^-Ui,.U^tJf^* 



a 



Airrif, 
^ = — 2 T 

01 02 



Im|f/f^^f/^^J/i* + tan/3A 



5Za = nib 



^{^Zl^^ + 5Zi^^ + 5Zaa) + sin /3 cos /3 5tan/3 



(78) 
(79) 

(80) 



5Z, 



terrif, 
2Mws^ 

U^mifi + tan/3A^) - f/^^^f/^ (/i* + tanM 



(5Z, 



2M2, s^ 



+ zlm {Ui^Ul^ (/. + tan/3A*)} 5Z?^^ + zim {t/^^.f/^^ (/i + tan/3A*)} (5Z^^;^ 

1 
~ 2 



f^6n f^L (^^ - /" tan /3) - f/^^^ Ul^ {Ab - /i* tan /3) 



5Z 



AG 
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The renormalization condition reads P|[TT| 



ReA(0, m? , m? ) + ReA(0, m? m? 



which corresponds to the two conditions 



62' b2' 



Re { ReA(0, m?^ , '^^J + ReA(0, ml^,ml^] 
Im <{ ReA(0, m?^ , m?^ ) + R;3A(0, m| , m?J 



0, 
. 



(82) 

(83) 
(84) 



The conditions (i)-(iv), are sufficient to fix all independent parameters and their respec- 
tive counterterms. As in Sect. I4.5[ relating Eqs. fl60|l and fl6T]) . one derives Eqs. (l62l) and 
which can also be written in the form 



Re {Al — fi tan /?} 6mh = —rrih Re SA^ — R,e6S + Im f/+ Im 6Yh 
Im {Al — fi tan (3} 6mh = +mi, Im 6A1, — Im 55' + Re U^ Im 6Yb 



^5) 
^6) 



with 



6S = -mfe(/i (5tan/3 + (5/itan/3) - Ul^Ui^^{5ml^ - 6mlJ - (Ret/+ - «Imt/_) Re(5n , (87) 

U± = Ur m ± f/ft U-* . (88) 

=■= Oil 622 "12 b21 ^ ' 

5m? and (5m? are fixed by Eqs. 042p and 04U|) . 



The above four equations ([H3D, dUD, (ISSj) and (15^ . can be solved for ReM;,, Im(5/lb, 
InK^y^ and 6mf,. Though, we still consider lieSAh and ImSAh as independent counterterms 
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we first calculate Iin 5Yi, and 5mi, in dependence of Re 5Ai, and Im 5Ai, for economically solv- 
ing the systems of equations. The solution for Im SYi^ and 5mb is 



Sruf, 
ImSYb 






with 



dr = 2 tan /3 (Im U+ Im 6S - Re U. Re 6S) 



+ 2 Re 



le 



26 Zo + ReA(0, ml^.mlj + ReA(0, m^^, m^ j + (5M + (5Zd[/, 



di = -2 tan (3 (Re [/+ Im 5^ + Im U. Re 5^) 



+ 2Im 



Mws^ 



le 



26Z^ + ReA(0, ml^^m^ + ReA(0, ml^.mQ] + 6M + SZ^U, 



Cr = +2 tan (3 [Im f/+ Im {Al - /i tan /?} - Re f/_ Re { A^ - /i tan /?}] + 2 Re ?7,„ , 

Cj = -2tan/3 [Ref/+Im{A^ -/itan/3} + Im[/_Re{yl^ -/itan/3}] + 2Imf/m , 

fr = -2lmUY , 

f, = 2t^nP{\UiJ-\U-,J) 
and 

^- = U-,^U~;jAl tan/3 + /i) - f/^^^^^^^(Atan/3 + /i*) . 
From the Eqs. fl85|) and fl86|) we immediately obtain ^Af, as 

Re 5 Ah = — f+ Im SYi, Im U+ - Re 55 - ^m,, Re {At - u tan /3}1 , 

Im 5Ab = — f- Im 5Yb Re U_ + Im 55 + 5mb Im {At - u tan/3}] . 
mfc 

Finally the Z factors in A have to be determined. The following condition is used 

R^S^Jm[)-ReS^Jmp = (z = l,2). 

This condition results in the following Z factors 



SZr 



ReS,Jm?)-ReS,Jm?) 



b2' 



2 

mf 

bi 



2 
b2 



1,2) 



(89) 
(90) 

(91) 

(92) 

(93) 
(94) 
(95) 
(96) 

(97) 

(98) 
(99) 

(100) 
(101) 



which guarantees the IR finiteness of the renormalized vertex A [TT 



Another subtlety has to be explained here: due to the fact that we have infrared divergent 
C-functions at pi = in A{pi = 0,p^,p^), we must deal with vanishing Gram-determinants. 
Therefore we follow Ref. [39] (and references therein) and replace the corresponding C- 
functions by well behaving linear combinations of 5-functions. Details can be found in the 
appendix. 
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4.7 Parameter definition 

The input parameters in the h/h sector have to correspond to the chosen renormahzation 
scheme. We start by defining the bottom quark mass, where the experimental input is the 
SM MS mass [lO], 

mp(mb) = 4.2 GeV . (102) 

The value of m^^{fiji) (at the renormahzation scale fin) is calculated from m^^{mi,) at the 
three loop level following the prescription given in Ref. 
An "on-shell" mass is derived from the MS mass via 



m?« = mrW l + ^^^^^|- + 21n^^ 



(103) 



The DR bottom quark mass is calculated iteratively frono 



we have used 










Sm^' = 


-Re {rrih 


[^Uml) + J:^{ml)] 


+ 


Smr = 


-Re {rUb 


KK) + Sf(m,^)] 


div 



m, - ^^^^ U04j 

with an accuracy of |1 — (m^^)'^"V("^fe'^)^"'~"'^^l < 10~^ reached in the nth step of the iteration. 
The bottom quark mass of a special renormahzation scheme is then obtained from 

mb = mp + (5mp - 6mi, . (105) 



[Sf(m^) + Sf(m^)]} , 

+ [Sf(m^) + SfV^)],J , (106) 

and 6mi) as given in Sects. I4.1H4.6I The quantity A;, [IHlIll] resums the (9((as tan/3)") and 
C((Q!i tan/?)") terms and is given by 

A, = ^^^tan/3M3>*/(m[,m|,m|) + ^^ tan/3 A,>* /(m|, m|, |/.n (107) 

with 

^, , , ab ln(6/a) + ac ln(a/c) + be ln(c/6) ,^^^. 

I{a, 6, c) = — ^ . 108) 

[a — c){c — b){b — a) 

Here at is defined in terms of the top Yukawa coupling yt{rnt) = y2'mt{mt)/v as at{mt) = 

y1{mt)/{A'K) with v = l/y^V2G^ = 246.218 GeV and mt{mt) ^ mt/{l - j^at{mt) + 
-^as{mt)). Ms is the soft SUSY-breaking parameter for the gluinos, with the gluino mass 
given as rrig := \Ms\. 



In case of complex Ab the replacement (1 + Ah) -^ |1 + Ab| should be performed [42] 



18 



5 Renormalization scheme analysis 

5.1 Calculation of loop diagrams 

In this section we give the relevant details about the calculation of the higher-order cor- 
rections to the decay channels (IT|2l) . Sample diagrams are shown in Figs. Ej |3l Not shown 
are the diagrams for real (hard or soft) photon and gluon radiation (which, however, can 
become numerically very important). They are obtained from the corresponding tree-level 
diagrams by attaching a photon (gluon) to the electrically (color) charged particles. The 
internal, in a generical way depicted particles in Figs. |2l |3] are labeled as follows: F can be 
a SM fermion, a chargino or neutralino or a gluino, S can be a sfermion or a Higgs boson, 
V can be a photon 7, a Z or W^ boson or a gluon g. Not shown are the diagrams with a 
gauge boson (Goldstone G'^)-Higgs selfenergy contribution on the external Higgs boson leg 
that can appear in the decay ^2 — > hiH^ . On the other hand, in our calculation, the wave 
function corrections for ^2 — )■ hiW^ vanish as all the external particle fields are renormalized 
on-shell. 

The diagrams and corresponding amplitudes have been obtained with the program Feyn- 
Arts [l5]. The further evaluation has been performed with FormCalc ^6]. As regularization 
scheme for the UV- divergences we have used constrained differential renormalization [17], 
which has been shown to be equivalent to dimensional reduction [18] at the one-loop level [16] . 
Thus the employed regularization preserves SUSY [l9l[50]. It was checked that all UV- 
divergences cancel in the final result. 

The IR-divergences from diagrams with an internal photon or gluon have to cancel with 
the ones from the corresponding real soft radiation. In the case of QED we have included 
the soft photon contribution following the description given in Ref. [27]. In the case of 
QCD we have modified this prescription by replacing the product of electric charges by 
the appropriate combination of color charges (linear combination of Ca and Cp times as)- 
More details will be given in Ref. [25]. Using the sbottom masses at the one- loop level, see 
Sect, m we found cancellation beyond one- loop order of the related IR and UV divergences 
for the decay ^2 — ^ hH^ , and a cancellation, as required, at the one-loop level for the decay 

t2 ^ m+% 

For completness we show here also the formulas that have been used to calculate the 
tree-level decay widths: 

r--(t2 ^ hH^) = ^' V ' (^ = 1' 2) , (109) 

|C(t2,6*,Vr)PA=^/2/ 2 2 ^2^ 



where A(a;, y,z) = {x — y — zY — iyz and the couplings C{a, b, c) can be found in the 
FeynArts model files [51]. The bottom- Yukawa couplings generically are enhanced with 
tan/3. 



^Using tree- level masses yields a cancellation of IR divergences beyond one- loop order also for ^2 ^- biW~^. 
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Figure 2: Generic Feynman diagrams for the decay ^2 -^ biH^ {i = 1,2). F can be a SM 
fermion, a chargino or neutralino or a gluino, S can be a sfermion or a Higgs boson, V can 
be a 7, Z, W^ or (7. Not shown are the diagrams with a W~^-H~^ or G^-H'^ transition 
contribution on the external Higgs boson leg. 



h 




s .^■ 







6,- 



t2 



s .^■ 



w 





t2 



I I 



V 



V 




w 



Figure 3: Generic Feynman diagrams for the decay t2 — ^ hW^ {i = 1,2). F can be a SM 
fermion, a chargino or neutralino or a gluino, 5* can be a sfermion or a Higgs boson, V can 
be a 7, Z, 1^=*= or g. 
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Seen. 


Mh± 


"^i2 


^^ 


At 


A, 


Ml 


Ma 


Ms 


SI 


150 


600 


200 


900 


400 


200 


300 


800 


S2 


180 


900 


300 


1800 


1600 


150 


200 


400 



Table 2: MSSM parameters for the initial numerical investigation; all parameters are in GeV. 
We always set rn^^imi,) = 4.2 GeV. In our analysis we use Mg^it) = Mi 



Ml 

Or 



=:M 



SUSY, 



where Msusy is chosen such that the above value of mi^ is realized. For the h sector the shift 
in Mq {b) as defined in Eq. (jH]) is taken into account. The parameters entering the scalar 
lepton sector and/or the first two generations do not play a relevant role in our analysis. The 
values for At and Ah are chosen such that charge- or color-breaking minima are avoided [52] . 



5.2 Numerical examples for the six renormalization schemes 

We start our analysis by showing some representative numerical examples. We evaluate 
the tree-level results and the one-loop correction for r(t2 — ^ biH'^) including wave function 
corrections. The parameters are chosen according to the two scenarios, SI and S2, shown in 
Tab. [210 

So far we concentrate on the rMSSM: if a scheme shows deficiencies in the rMSSM, 
the same problems occur in the cMSSM. The final numerical examples in Sect. [6] will also 
show complex parameters as well as results for r(t2 — ^ &i,2W^^)- It should be noted that 
tan/3 < 9.6 (4.6) is excluded for SI (S2) due to the MSSM Higgs boson searches at LEP [531 
^^. However, we are interested in the general behavior of the renormalization schemes. If 
certain features appear in the two numerical scenarios (SI and S2) only for experimentally 
excluded tan /3 values, other parameter choices may exhibit these features also in unexcluded 
parts of the MSSM parameter space. Consequently, in order to investigate the various 
renormalization schemes on general grounds, in the following we show the results for tan /3 > 
1. A similar reasoning applies to the limits on the MSSM parameter space due to SUSY 
searches. Nevertheless, to avoid completely unrealistic spectra, the following exclusion limits 
hold in our two scenarios: 



m 



;^ > 95 GeV, m^^ > 89 GeV, niq > 379 GeV, rrie, > 73 GeV, 



m^o > 46 GeV, m-± > 94 GeV, m~„ > 308 GeV. 



"Xi 



fiir 



A few examples of the scalar top and bottom quark masses at the one-loop levejj (using 
M? (6) in Eq. (jH]) for the one-loop result) in the scenarios SI and S2 are shown in Tab. [31 



^It should be noted that we do not include any further shifts in the parameters than the one given in 
Eq. p4)) . Correspondingly, the values for the parameters Ai, and Afg in Tab. [2] do not reflect the actual 

values for the input parameters with respect to the chosen renormalization scheme. For example, the 62 mass 
— though considered as an input in the renormalization scheme and defined as on-shell mass — receives a 
shift going from tree- to one-loop level when starting out with the values in Tab. [2] and including only the 
shift Eq. (j44| . To circumvent this shift of the 62 mass, additional shifts to the tree-level values of Ab and 
Mg would be required (depending on the renormalization scheme) . 

^For the scalar top quark masses the tree-level and the one-loop values are the same (according to our 
renormalization conditions) . 
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Seen. 


tan/3 


m. 


m^ 


^l. 


^l. 


SI 


2 


293.391 


600.000 


441.987 


447.168 


20 


235.073 


600.000 


418.824 


439.226 


50 


230.662 


600.000 


400.815 


449.638 


S2 


2 


495.014 


900.000 


702.522 


707.598 


20 


445.885 


900.000 


678.531 


695.180 


50 


442.416 


900.000 


628.615 


697.202 



Table 3: The top and bottom squark masses at the one-loop level (see text) in the scenarios 
SI and S2 and at different tan/3 for the numerical investigation; all masses are in GeV and 
rounded to one MeV. 



The values of mi^ allow copious production of the heavier scalar top quark at the LHC. For 
other choices of the gluino mass, rrig > mi^ , which would leave no visible effect for most of 
the decay modes of the ^2, the heavier scalar top quark could also be produced from gluino 
decays at the LHC. Furthermore, in SI (even for the nominal value of mi^ as given in Tab. [2]) 
the production of ^2 at the ILC(IOOO), i.e. with y/s = 1000 GeV, via e+e" — )■ t2^i will be 
possible, with the subsequent decay modes ([I]) and ([2]) being open. The clean environment 
of the ILC would permit a detailed study of the scalar top quark decays. Depending on the 
combination of allowed decay channels a determination of the branching ratios at the few 
per-cent level might be achievable in the high-luminosity running of the ILC(IOOO). More 
details will be discussed elsewhere [25] . 

Later we will also analyze numerical results for complex input parameters. Here it should 
be noted that the results for physical observables are affected only by certain combinations of 
the complex phases of the parameters fj,, the trilinear couplings Af, f = {u, c, t, d, s, b, e, fi, r}, 
the gaugino mass parameters Mi, M2, M3 and the Higgs soft SUSY breaking parameter 
"^12 [55l[56]. It is possible, for instance, to eliminate the phase ipM2 a^id the phase ^9^2 . 
Experimental constraints on the (combinations of) complex phases arise in particular from 
their contributions to electric dipole moments of heavy quarks [57], of the electron and 
the neutron (see Refs. [58l|59] and references therein), and of the deuteron [60]. While SM 
contributions enter only at the three-loop level, due to its complex phases the MSSM can 
contribute already at one-loop order. Large phases in the first two generations of sfermions 
can only be accommodated if these generations are assumed to be very heavy [61] or large 
cancellations occur [62], see however the discussion in Ref. [631 IM] . A recent review can be 
found in Ref. [65]. Accordingly, using the convention that ipM2 = and ^9^2 = 0, as done 
in this paper, in particular the phase ip^ is tightly constrained [61], while the bounds on the 
phases of the third generation trilinear couplings are much weaker. The phase of /i enters 
in the combinations {(pAtb + V^^ ~ '^m'^ )• Setting y?^ = (and (pM2 = ^m^ = 0, see above) 
leaves us with At and Ab as complex valued parameters. Since we are interested in the 
renormalization of the b/b sector, in our numerical analysis we will focus on a complex Ai, 
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r(t2 ^ biH+) for SI 


tan/3 = 2 


tan /3 = 50 


renorm. scheme 


^^R 


tree 


loop 


rrib 


tree 


loop 


rrib 


RSI: "OS" 


m. 


0.0017 


-0.0011 


3.29 


2.5930 


-53.3469 


3.84 


RS2: "mfe, A^ DR" 


m^ 


0.0009 


0.0002 


2.38 


0.9653 


-0.0311 


2.16 


RS3: "mb, Yt DR" 


m. 


0.0009 


0.0004 


2.38 


0.9484 


-1.5404 


2.16 


RS4: "mft DR, n OS" 


^h 


0.0009 


0.0000 


2.38 


0.9593 


-0.3411 


2.16 


RS5: "ADR, Ren OS" 


m. 








0.9399 


-0.0481 


2.13 


RS6: "A vertex, Ren OS" 


^h 


0.0007 


0.0001 


2.19 


0.9390 


-0.0347 


2.13 




r(t2 -^ hiH+) for S2 


tan/3 = 2 


tan /3 = 50 


renorm. scheme 


l^R 


tree 


loop 


rrib 


tree 


loop 


rrib 


RSI: "OS" 


^2 


2.0928 


-0.0776 


3.23 


8.5163 


-106.9700 


3.70 


RS2: "mft. At, DR" 


^h 


2.2171 


-0.1449 


2.33 


1.8173 


-0.5125 


2.11 


RS3: "mb, n DR" 


m. 


0.0077 


0.0582 


2.33 


3.1409 


-11.6833 


2.11 


RS4: "mfe DR, n OS" 


^h 


2.2564 


-0.1031 


2.33 


2.9230 


-4.5506 


2.11 


RS5: "ADR, Ren OS" 


^h 


2.2332 


-0.1004 


2.45 


2.3018 


0.2924 


1.84 


RS6: "A vertex, Ren OS" 


^h 


2.2925 


-0.1067 


2.14 


2.3558 


-0.0710 


1.86 



Table 4: Examples for tree-level and full one-loop contributions (see text) to r(t2 — >■ biH"^) 
for SI (upper table) and S2 (lower table); all values are in GeV (no comparison of the 
renormalization schemes, see text). In SI using RS5 a divergence is reached for tan/3 = 
|/lfe|/|/i| = 2 and no value can be computed (see text below). The different renormalization 
schemes are listed in Tab. [1] 



and keep At real (see, however, Ref. [25]). 

We start our numerical examples with the evaluation of r(t2 — ^ biH~^) in SI and S2 for 
tan/3 = 2 and tan/3 = 50 as shown in Tab. HI The corresponding results as a continuous 
function of tan /3 can be seen in Fig. |H It must be emphasized here that the table and 
the plots do not constitute a comparison of the various schemes, but "only" individual 
numerical examples that are used to exhibit certain problems of the various schemes. A 
numerical comparison of the schemes requires that the input parameters are converted from 
one scheme into another, see, for instance, Ref. [10], which is not performed within this 
analysis. In our numerical examples the renormalization scale, fiR, has been set to the 
mass of the decaying particle, i.e. ^r = mi^. In Tab. H] the two main columns, labeled 
"tan/3 = 2" and "tan/3 = 50", are divided into three columns where "tree" contains the 
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tan/3 



10 15 20 25 30 35 40 45 50 



10 15 20 25 30 35 40 45 50 



Figure 4: r(t2 -^ biH'^). Full one-loop corrected partial decay widths for the different 
renormalization schemes (no comparison, see text). The parameters are chosen according to 
SI in the left plot and S2 in the right plot. For SI the grey region and for S2 the dark grey 
region is excluded by LEP Higgs searches (see text). 



tree-level results and "loop" the one-loop contribution, rrif, denotes the corrected bottom 
quark value corresponding to the respective renormalization, see Eq. (llOSp . 

The two values of tan/3 were chosen as an example of a very low and a very high value. 
It should be kept in mind that the low value is possibly already in conflict with MSSM 
Higgs boson searches [S31[M], but kept to show an "extreme" example as explained above. 
It can be seen that RSI, RS3, RS4 and RS5 yield relatively large absolute values of loop 
contributions with respect to the tree-level result, either for tan/3 = 2 or for tan/3 = 50, at 
least in one of the two numerical scenarios. This simple example shows that (by choosing a 
specific scenario) already all except two renormalization schemes fail in part of the parameter 
space. 

More problems of the renormalization schemes RSI, RS3, RS4 and RS5 become visible in 
Fig. m In the left (right) plot of Fig. |4]we show the results of S1(S2) as a function of tan/3. 
For SI the grey region and for S2 the dark grey region at low values of tan/3 are excluded by 
LEP Higgs searches [51]. It can be seen in Fig. |4]that RSI and RS3 deviate strongly from 
the (see the end of Sect. 15. ip expected behavior of increasing r(t2 — ^ biH'^) with growing 
tan/3 that the other schemes exhibit. The same is observed for RS4 in S2 for tan/3 > 35. 
Problems in RS2 are discussed in Sect. 15. 6^ problems in RS6 have been found for complex 
parameters, see Sect. 15.71 The various spikes and dips can be understood as follows: 



For RS3 in S2 a "peak" appears at tan/3 
in Sect. 15.51 below. 



4.6 and at tan/3 ^ 6.2. This is discussed 



For RS5 in SI a "peak" appears (not visible) at tan/3 = |ylh|/|/i| = 2. This is caused 
by large corrections to the bottom quark mass as discussed further in Sect. 15.71 This 
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is also the reason why there is no entry in Tab. IHfor RS5, SI at tan/3 = 2. 

• For RS5 in S2 a "peak" appears at tan/? = |Ab|/|yu| = 5.33. This is caused by large 
corrections to the bottom quark mass as discussed further in Sect. 15.71 

5.3 Generic considerations for the 6/6 sector renormalization (I) 

As discussed in Sect. HI a bottom quark/squark sector renormalization scheme always con- 
tains dependent counterterms which can be expressed by the independent ones. According 
to our six definitions, this can be 6mh, 6Ai, or 6Yi,. A problem can occur when the MSSM 
parameters are chosen such that the independent counterterms (nearly) drop out of the re- 
lation determining the dependent counterterms. As will be shown below, even restricting 
to the two numerical examples, SI and S2, it is possible to find a set of MSSM parameters 
which show this behaviour for each of the chosen renormalization schemes. Consequently, 
it appears to be difficult by construction to define a renormalization scheme for the bottom 
quark/squark sector (once the top quark/squark sector has been defined) that behaves well 
for the full MSSM parameter space. One possible exception could be a pure DR scheme, 
which, however, is not well suited for processes with external top squarks and/or bottom 
squarks. 

Assuming that SUSY, and more specifically the MSSM, will be discovered at the LHC 
and its parameters will be measured, the problem will have disappeared. For a specific 
set of MSSM parameters, renormalization schemes can (easily) be found that behave well. 
However, due to our ignorance about the actual values of the SUSY parameters, scans over 
large parts of the MSSM parameter space are performed, see also Sect. El For this kind of 
analysis a careful choice of the renormalization scheme has to be made. 

In the following subsections we will analyze in more detail, analytically and numerically, 
the deficiencies of the various schemes. 

5.4 Problems of the "OS" renormaUzation 

The "OS" renormalization as described in Sect. 14.11 does not yield reasonable results in 
perturbative calculations as shown already, e.g., in Ref. [TOl[TT]. For the sake of completeness 
we briefiy repeat the results. The "OS" scheme of Sect. 14.11 is the renormalization scheme 
analogous to the one used in the t/i sector and thus would be the "naive" choice. It includes 
an on-shell renormalization condition on the sbottom mixing parameter Yf, that contains the 
combination [Ah — yU*tan/3). In parameter regions where (/itan/3) is much larger than Af,, 
the counterterm 6Ab receives a very large finite shift when calculated from the counterterm 
6Yb. More specifically, 6Ab as given in Eq. ( 148|) contains the contribution 

dAb = — [-(Ab- fx* tan (3) Smb + ...] (112) 

nib 

that can give rise to very large corrections to Ab. This is also visible in Fig. [5]below, where we 
show the numerical values of 6Ab as a function of tan /3 for various renormalization schemes. 
In Ref. [To] it was shown that, because of Eq. (11121) . the "OS" renormalization yields huge 
corrections to the lightest MSSM Higgs mass. Also the numerical results shown in Tab. H] 
and Fig. |3]show extremely large one-loop corrections for tan/3 = 50. 
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Finite parts of SAi, in various renormalization schemes. The parameters are chosen 
to SI left plot and S2 right plot. For SI the grey region is excluded and for S2 the 
region is excluded. 



This problem is (more or less) avoided in the other renormalization schemes introduced in 
Tab. [T], where the renormalization condition is applied directly to Af,, rather than deriving 
SAf) from a renormalization condition fixing SYf,. Also the renormalization schemes RS3 
("mb, Yb DR") and RS4 ("rrab DR, Yh OS") avoid this severe problem by renormalizing the 
bottom quark mass DR. 

5.5 Problems of non-Ab renormalization schemes 

Two of our schemes, besides the "OS" scheme (RSI), do not employ a renormalization of Ab. 
RS3 ("mfc, n DR") and RS4 ("mj, DR, ^ OS"). As argued in Sect. [El a huge contribution 
to 6Ah as evaluated in that section is avoided by the DR renormalization of rrib. However, 
following Eq. fHHl) with 5m? , 5m? , 6Yi, and 6mi, chosen according to the renormalization 
schemes RS3 and RS4, respectively, one finds for the finite parts of 6Ai,: 



RS3: 6 A, 



fe fin 



RS4: MJfin 



nib L 
1 

nib 



Ur m (5nil -5ml 

Oil 6i2 y bi 62 

Ur W (6ml -6ml 

Oil 6l2 V 01 62 



fin 



Ui W 6Y: + UI Ur 6Yb 

Oil b22 " fci2 021 " 



fin 



+ 



(113) 
(114) 



where the ellipses denote contributions from 6fi which, however, are not relevant for our 
argument. It can be seen that still 6Ab depends on parameters (diagonal and off-diagonal 
sbottom self-energies) that are independent of Ab. As an example, Higgs boson loops in the 
sbottom self-energy contain contributions ~ /itan/3, which can become very large, indepen- 
dently of the value of Ab. This can be seen in the right plot of Fig. [5], where we show 6Ab as 
a function of tan/3 in S2. In both renormalization schemes, RS3 and RS4, 6Ab becomes very 
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large and negative for large tan /3. This yields the very large and negative loop corrections 
to r(t2 -^ biH^) shown in the right plot of Fig. HI In SI this problem is less pronounced, as 
can be seen in the left plot of Fig. [5] (SAb) and Fig. H (1(^2 -^ biH~^)). 

But also for lower tan/3 values, tan/3 ^ 10, problems can occur. The (finite) "multiple 
spike structure" in RS3 for S2 around tan/3 f» 5.33 (for details see the small insert within the 
right plot of Fig. |5]) is due to an interplay of top/chargino contributions to the two diagonal 
sbottom self-energies, invalidating this scenario also for this part of the parameter space. 

5.6 Problems of an mb—Ab renormalization 

If rrib and A^ are renormalized, the sbottom mixing parameter Yf, is necessarily a dependent 
parameter, see Tab. [1] This situation is realized in the scheme RS2 ("m^, A^ DR"), see 
Sect. 14.21 6Yh enters prominently into SZi^_^. For real parameters we have. 

Re Ez {ml ) - 6Yb 
SZi,, = -2 '-^^^ . (115) 

21 j^Z _ ^ 

01 02 

In this way SY^ (or the interplay between SYh and ReSj2^(m? )) can induce large loop cor- 
rections to the scalar top quark decay width. 6Yh can be decomposed according to Eq. (ICTl 
(concentrating again on the case of real parameters), 

^^^ = i^,'[^'t^, p (^-ii - ^<) + • ■ ■ ' (116) 

l^felll I^bl2l ^ ^ 

where the ellipses denote terms with only divergent contributions (due to the chosen renor- 
malization scheme RS2) as well as finite contributions from 6fi, which, however, do not play 
a role for our argument. For "maximal sbottom mixing", |f/^ | si |[/^ |, SYb diverges, and 
the loop calculation does not yield a reliable result. In our two parameter scenarios, SI and 
S2, this is not the case. Such a large sbottom mixing is often associated with large values of 
\Ab\ that may be in conflict with charge- or color-breaking minima |52j . 

However, in order to show an example with a divergence in 6Yb we use a modified version 
of SI with Ab = 1000 GeV (a value still allowed following Ref. [S2]). In this scenario at 
tan/3 ^ 37 we indeed find the case of "maximal mixing" in the scalar bottom sector. As 
expected this leads to a divergence in SYb, as can be seen in the left plot of Fig. [61 This 
divergence propagates into 6Zi^_^ as shown in the right plot of Fig. [6]cl(Also E^^^ exhibits a 
discontinuity due to a sign change in Ui for this extreme set of MSSM parameters.) The 
tan/3 value for which this "divergence" occurs depends on the choice of the other MSSM 
parameters. For (numerical) comparison we also show SZi^_^ for the two scenarios. 

For the different choice of MSSM parameters in S2 (without a higher Ab value) this 
divergences does not occur. However, for tan /3 ^ 7 one finds 6Yb ^ Re S^^^ {m? ) (with 
SYb = ReS^ (m? ) for tan/3 ^ 7.5). In this part of the parameter space we also find 
m^^ ~ TTi^^, yielding a relatively large value of ^^^^i according to Eq. (11151) . as can be seen 
in the right plot of Fig. O This relatively large (negative) value of SZi^_^ in turn induces 



^ The scalar bottom masses could receive large corrections via Af ? (&) in Eq. pi)) . with SYi, entering via 
Eq. (gSl). 
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Figure 6: T(t2 — >■ hiH^). Left plot: size of SY}, and ReS^^^(m? ), the two contributions to 
5Z^^^, in RS2. Right plot: comparison of the size of ^Z^^^ in the scheme RS2 ("m;,, A;, DR"). 
For both plots the parameters are chosen according to SI (but here with Ab = 1000 GeV), 
S2 in Tab. [2l For SI the grey region is excluded and for S2 the dark grey region is excluded 
via LEP Higgs searches (see text). 



relatively large corrections to T{t2 — )■ hiH^). However, the loop corrections do not exceed 
the tree-level value of r(t2 — ^ biH~^) (for our choice of MSSM parameters). In summary: 
while for SI a divergence in 6Yh and thus in SZ^ can appear for very large values of \Ah\ 
(possibly in conflict with charge- or color-breaking minima), invalidating the renormalization 
scheme RS2 in this part of the parameter space, these kind of problems are not encountered 
in S2. Here only moderate loop corrections to the respective tree-level values are found, and 
RS2 can be applied safely. 



5.7 Problems of non-mb renormalization schemes 



Two of our schemes do not employ a renormalization condition for nif,: RS5 ("A;, DR, 
ReYj, OS") and RS6 ("Afe vertex, ReF^ OS"). Since A^ and Yj, are complex, we chose to 
renormalize A), and the real part of Y^. 

We start with the discussion of the (simpler) "Af, DR, Re Yh OS" scheme. We will focus on 
the real case as a subclass of the more general complex case. In this renormalization scheme 
the bottom quark mass counterterm has the following form for real parameters (compare to 
Eq. dSSD), 



drrih 



ruh SAh + SS 
{Ab — /itan/3) 



:il7) 



For vanishing sbottom mixing one finds {Ab—fi tan /3) — ?■ 0. In the "Ah DR, Re Yb OS" scheme 
this yields a finite (and negative) numerator in Eq. flll7p . but a vanishing denominator. 
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Figure 7: Left plot: m^ in RS5 and RS6 for SI, S2. For SI the grey region is excluded and 
for S2 the dark grey region is excluded. Right plot: m^ in RS6 for SI, S2 but both with 
tan/3 = 20 and ipA,, varied. In S2 we used also |/i| = 120 GeV. 



In a numerical evaluation, starting out with a value for the bottom quark mass defined 
as DR parameter, the actual value of the bottom quark mass receives a shift with respect to 
the DR bottom quark mass according to Eq. f llOSp . This shift corresponds to the finite part 
of Srrib in Eq. flll7p . Consequently, large positive or negative contributions to the bottom 
quark mass can occur, yielding possibly negative values for the bottom quark mass and thus 
invalidating the renormalization scheme for these parts of the parameter space. This can be 
seen in the left plot of Fig. [71 where we show rrih in RS5 (and RS6) for the two numerical 
scenarios given in Tab. [2] as a function of tan/?, rrib exhibits a strong upward/downward 
shift around the pole reached for tan/3 = A^/n and consequently yields unreliable results in 
this part of the parameter space. 

We now turn to the RS6 ("Ab vertex, lie Yh OS") scheme. Following the same analysis 
as for the "A;, DR, Re Yj, OS" scheme an additional term in the denominator of the bottom 
quark mass counterterm ~ Um/U^ appears, 

dS + F 



Sirih 
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(Afe-/itan/3) - Um/ {tan PU-) 

where F denotes other (relatively small) additional contributions. With the help of Eq. (|9 
one finds for real parameters 

Um _ Ui^, Ui^, {A tan (3 + fi)- U-,^^ U-,^^ {A^ tan /3 + /i) 



tan/3[/_ 
and therefore 



tan P{U,^U,^^-U,^U-, J 

6S + F 
onih 



fi (tan/3 + l/tan/3) 



(A + /i/tan/3) , (119) 



(120) 
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The denominator of Eq. (11201) can go to zero only for fi ^ 0, which is experimentally already 
excluded. Consequently, the problem of (too) large contributions to rrih is avoided in this 
scheme. This can be seen in the left plot of Fig. [71 where RS6, contrary to RS5, does not 
exhibit any pole-like structure in rrib. 

In the complex case the above argument is no longer valid, and larger contributions to 
Srub can arise. In the limit of tan /3 ^ 1 and n real the denominator of 6mh in Eq. fl89|l reads 

-^~4/itan3/3[Re[/_(|f/r \^ - \Ur \^) +lmU_ .^"^^ „ Im fa K a) 1 . (121) 

Depending on 99^^ this denominator can go to zero and thus yield unphysically large correc- 
tions to m^ in RS6. In the right plot of Fig. [7]we show TTi^ as function of ifA^- At y?^^ ~ 41.5°, 
87.5°, 272.5°, 318.2° the denominator in Eq. (I12ip goes to zero and changes its sign which 
explains the corresponding structures. This divergence in 5mi, enters via Eq. (I105p already 
into the tree- level prediction. To summarize: while in SI the scheme RS6 is well-behaved 
and can be safely applied (also for complex Ab), in S2 (with |/i| = 120 GeV) severe problems 
(divergences in the counterterms) arise once complex parameters are taken into account. 
Consequently, for S2 the scheme RS6 cannot be applied. 

It should be noted that the "Af, vertex. Re Y^ OS" (RS6) scheme is the complex version of 
the renormalization scheme used in Refs. [IDIIII] for the 0{abOts) corrections to the neutral 
Higgs boson self-energies and thus to the mass of the lightest MSSM Higgs boson, M^. For 
real parameters, no problems occured. Therefore, employing this renormalization scheme in 
Refs. [TOllTT] yields numerically stable results. 

5.8 Generic considerations for the 6/6 sector renormalization (II) 

In the previous subsections we have analyzed analytically (and numerically) the deficiences of 
the various renormalization schemes. We have shown that despite of the variety of schemes, 
even concentrating on the two sets of parameters, SI and S2, severe problems can be en- 
countered in all schemes. 

For the further numerical evaluation of the partial stop quark decay widths we choose 
RS2 as our "preferred scheme". According to our analyses in the previous subsections, 
RS2 shows the "relatively most stable" behavior, problems only occur for maximal sbottom 
mixing, \Ui \ = \Ui |, where a divergence in 6Yf, appears. Having 6Yb as a dependent 
counterterm induces effects in the field renormalization constants SZ^ and SZ^ and in 
5M? (b) entering the scalar bottom quark masses. In a process with only internal scalar 
bottom quarks, no problems occur due to the field renormalization, but counterterms to 
propagators, which induce a transition from a bi squark to a 62 squark contain also the 
term 6Yb. However, 6Yb appearing in counterterms of internal scalar bottom quarks does 
not exhibit a problem, since in this case these "dangerous" contributions cancel (which we 
have checked analytically). On the other hand, other schemes with 6m,b or 6Ab as dependent 
counterterms may exhibit problems in larger parts of the parameter space and may induce 
large effects, since m,b (or the bottom Yukawa coupling) and Ab enter prominently into the 
various couplings of the Higgs bosons to other particles. 

We are not aware of any paper dealing with scalar quark decays (or decays into scalar 
quarks) that has employed exactly RS2 (or its real version), see our discussion in the be- 
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ginning of Sect. HI Very recently a calculation of the scalar top decay width in the rMSSM 
using a pure DR scheme for all parameters was reported [66] . 



6 Numerical examples for our favorite scheme 

Following the discussion in Sect. [5] we pick the renormalization scheme that shows the "most 
stable" behavior over the MSSM parameter space. We choose the "mfe, Af, DR" (RS2) scheme. 
Tree-level values of the partial decay widths shown in this section have been obtained in- 
cluding a shift in rrih according to Eq. (11051) . We will concentrate on the calculation of the 
partial ^2 decay widths including one scalar bottom quark in the final state. A calculation of 
the respective branching ratios requires the evaluation of all partial scalar top quark decay 
widths, which in turn requires the renormalization of the full cMSSM. This is beyond the 
scope of our paper and will be presented elsewhere fI5\ . 

6.1 Full one- loop results 

We start our numerical analysis with the upper left plot of Fig. [HI where we show the partial 
decay width T(t2 -^ hiH^) as a function of tan/3, "tree" denotes the tree-level value and 
"full" is the decay width including a// one-loop corrections as described in Sect. 15.11 As one 
can see, the full one-loop corrections are negative and rather small over the full range of 
tan/3, the largest size of the loop corrections is found to be ~ 28% of the tree-level value for 
tan/3 = 50 in S2|3 In SI the grey region and in S2 the dark grey region is excluded due to 
too small values of the mass of the lightest MSSM Higgs boson, Mh- 

In the upper right plot of Fig. [8] we show the partial decay width varying \Ai,\ for tan/3 = 
20. In SI and S2 the full one-loop corrections grow with A?,, but never exceed ~ 25% of the 
tree-level result. Note, that for SI \A^\ > 1130 GeV (grey region) and S2 \Ab\ > 1800 GeV 
(dark grey region) is excluded due to the charge- or color-breaking minima. Over the full 
parameter space the loop corrections are smooth and small with respect to the tree-level 
results. 

In the lower left plot of Fig. [S]we analyze the partial decay width varying |/i| for tan/3 = 
20. Values for |yu| < 120 GeV are excluded due to m^± < 94 GeV [IQ]. The loop corrected 
predictions for the partial decay width show several dips and spikes. In SI the first dip at 
|/i| ^ 285 GeV is due to If/^^J ~ l^bijl' ^^^ ^^^ discussion in Sect. 15.61 The second peak/dip 
(already present in the tree- level prediction) at |/i| = 300 GeV is due to the renormalization 
of fi [12] and will be discussed in more detail in Ref. [23] jj The third dip at |/i| ^ 424 GeV, 
which is hardly visible, is due to the production threshold rrit + m^o = mi^. The fourth 
dip at |/i| ~ 873 GeV is the threshold m^-^ -|- Mh± = rni^ of the self energy S^^^(m? ) in the 
renormalization constants dZ^ and drrtj^ . The fifth dip at |/i| ~ 1107 GeV is the production 
threshold m^^ + M^ = m^^. For |/i| > 790 GeV the value of M^ drops strongly, and the 
scenario SI is excluded by LEP Higgs searches as indicated by the gray shading. Apart from 



"^ It is interesting to note that at tan/3 = |Af,|/|/.i| = 2 (5.33) in SI (S2) we get t/g =1 and t/g =0, 

and consequently bL,R = &i,2, respectively. 

^ The chosen renormalization exhibits a divergence for /i = M2. 6fi enters via SYi, into SM^ (b) and thus 
into the values of mj^ . Consequently, the dip is already present in the tree-level result. 
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the dips analyzed above the loop corrections are very small and do not exceed ~ 7% of the 
tree- level result, the prediction for T{t2 — t- biH'^) is well under control. We now turn to the 
scenario S2. Here, for growing |/i|, the squark mass splitting in the t/b sector becomes very 
large, leading to large contributions to the electroweak precision observables. The dark gray 
region for |/i| > 1060 GeV yields W boson masses outside the experimentally favored region 
at the 2(7 level, Myy ^ 80.445 GeV [67]. Such large |/i| values are consequently disfavored. 
The dip/peak at |/i| = 200 GeV in the tree and the loop contribution is due to 6fi, where 
H = M2 is reached, see above. The second dip at |/i| = 477 GeV, which is hardly visible, is 
the threshold nit + ''^x^ = ''^bi °^ ^^^ ^^^^ energy S^^^ (?7i? ) in the renormalization constants 
SZr and 6ml . The third dip at |/i| = 725 GeV is the production threshold mt + m^o 
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mz. 



The fourth dip at |yu| = 850 GeV is again the threshold tti^-^ +Mh± = t^i^- In S2 the one-loop 
corrections are negative and growing with |/i|. Apart from the dips described above, also in 
this numerical evaluation the loop corrections stay mostly relatively small with respect to 
the tree- level result, reaching the largest relative contribution at the smallest |/i| values, and 
are thus well under control. 

We now turn to the case of complex parameters. As discussed in Sect. 15.21 we consider 
only Ab as a complex parameter. In the lower right plot of Fig. [8] we show the partial decay 
width depending on Lp^^ for tan/3 = 20. In SI, the tree-level values and the loop corrections 
are well-behaved. The latter ones stay relatively small for the whole parameter space, not 
exceeding ~ 18% of the tree-level result. In S2, the largest corrections occur for real positive 
values of A^ and reach ~ 12% of the tree-level values. For negative A^, the tree-level result 
becomes very small (< 0.01 GeV) and here the size of the loop corrections can be as large 
as the tree-level values. A small (and barely visible) asymmetry in the one-loop corrections 
appears in the lower right plot of Fig. [HI due to terms ~ f/r x Co,i,2-function. The peak/dip 
at ifA^ ~ 117°, 243° are again due to If/^nl ~ 1^612!' ^^^ Sect. 15.61 It can be seen that the 
peaks due to this divergence are relatively sharp, i.e. the region of parameter space that is 
invalidated remains relatively small. 

In Fig. [9] we show the results for r(t2 — ^ &2-f^^) for the same set and variation of pa- 
rameters as above. Consequently, the same peak and dip structures are visible in the lower 
plots of Fig. [91 In the lower left plot of Fig. [9] in SI both lines end because the phase space 
closes, mi^ +Mh± > mi^ for |yu| > 300 GeV. Overall the partial decay width is much smaller 
than for r(t2 — >■ biH~^), which can partially be attributed to the smaller phase space, see for 
instance the results within S2 in the upper left plot of Fig. [9], and partially to the smallness 
of the tree-level coupling. Only in S2 for tan/3 ^ 35 we find r(t2 — )■ &2-f^^) ~ 1 GeV. The 
relative corrections become very large for |v4fe| ^ 1200 GeV as shown in the upper right plot 
of Fig. [9], however these values are disfavored by the constraints from charge and color break- 
ing minima as discussed above. The smallness of r(t2 — ^ &2-f^^) at the tree-level can lead 
sometimes to a "negative value at the loop level". In this case of (accidental) smallness of 
the tree- level partial decay width also |A^ioopP would have to be taken into account, yielding 
a positive value for T{t2 — ?■ b2H^). Overall, because of the smallness of the tree-level result 
due to the tree-level coupling the relative size of the loop corrections are a bit larger than 
for r(t2 — ^ biH^). Nevertheless, apart from the peaks visible in the lower plots of Figs. [9], 
the loop corrections are well under control also for r(t2 — ^ b2H^) using the renormalization 
scheme RS2. Again a small asymmetry in the one-loop corrections in the lower right plot of 
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Fig. [9] can be observed, which is due to terms ^ Ur,, x Cq i 2-function. 

Finally we evaluate the partial decay width of a scalar top quark to a scalar bottom quark 
and a W boson, r(t2 — > biW~^) and r(t2 — > b2W^). Since the W boson is relatively light, 
also the latter channel is open. In Fig. [10] the results for ^2 ~^ biW^ are shown, in Fig. [11] 
the ones for ^2 ~^ b2W^. The divergences visible in the various plots are the same ones as 
found in the respective plot for r(t2 — ^ biH^). An additional (finite) dip is visible in the 
lower left plot of Fig. [TO] in S2 for |/i| '^ 521 GeV, due to an interplay of t/x2 contributions 
to S^^^ ("^bj, similar to the structure discussed for Fig. [5] In this part of the parameter space 
the results calculated within the renormalization scheme RS2 have to be discarded. 

Overall, the loop corrections to r(t2 — ^ biW^) calculated within the renormalization 
scheme RS2 behave similar to the ones to r(t2 — ^ biH^). The size is relatively small, i.e. 
^ 20% and ^ 30% of the tree-level results in the upper left and in the upper right plot of 
Fig. [ini respectively, for the regions which are not in conflict with charge- or color breaking 
minima (for |Afe| = 2000 GeV a correction of ~ 70% of the tree-level result can be observed 
in SI due to the smallness of the tree-level value). We find loop corrections of the size of 
^ 20% of the tree-level results in the lower left plot of Fig. [10] except for very small values 
of |/x| and in the lower right plot of Fig. [TU] In the latter plot for S2 the known divergences 
appear at lpa^, ~ 117°, 243°, leading to larger loop corrections for intermediate values of ipA,,- 
Apart from the latter case the full one-loop corrections to r(t2 — ^ biW^) are well under 
control employing the renormalization scheme RS2. 

Similar observations hold for the decay ^2 -^ b2W~^, as shown in Fig. [11] In the upper left 
plot of Fig. [Uin the scenario S2 for tan/3 = |74ft|/|/x| (^ 5.3, the tree-level partial decay width 
vanishes, leading to a "negative value at the loop level" . As discussed above, in this case also 
lA^ioopP would have to be taken into account, yielding a positive value for r(t2 -> b2W^). (A 
similar situation is found in the lower left plot of Fig. [11] for \^\ ^ 200 GeV.) For somewhat 
larger tan/3 values, loop corrections of ~ 50% of the tree-level values are reached, while in 
SI they stay below ~ 23% of the tree-level results. In the upper right plot of Fig. [TT] the 
loop corrections are smaller than ~ 40% of the tree-level values, depending on the size of 
|A(,|, see above. The loop corrections shown in the lower left plot of Fig. [TT] yield maximal 
~ 9(37)% of the tree-level results in SI (S2), apart from very small fi values, where the 
tree-level partial decay width can become accidentally small. 

Finally, looking at the dependence on ipAt in the lower right plot of Fig. [Til apart from 
the known divergences in S2 around (pA^ ~ 117°, 243°, the loop corrections do not exceed 
~ 6% and ~ 35% of the tree-level values in SI and in S2, respectively. Overall, except 
for the small parameter regions around ipA^ ~ 117°, 243°, the full one- loop corrections to 
r(t2 — ^ 62 W^^) are well under control employing the renormalization scheme RS2. 
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Figure 8: r(t2 — ?■ biH^). Tree-level and full one-loop corrected partial decay widths for 
the renormalization scheme RS2. The parameters are chosen according to the scenarios SI 
and S2 (see Tab. |2]). For SI the grey region is excluded and for S2 the dark grey region is 
excluded. Upper left plot: tan/3 varied. Upper right plot: tan/3 = 20 and \Ai,\ varied. Lower 
left plot: tan/3 = 20 and |/x| varied. Lower right plot: tan/3 = 20 and ipA^ varied. 
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Figure 9: r(t2 — ^ &2-f^^)- Tree-level and full one-loop corrected partial decay widths for 
the renormalization scheme RS2. The parameters are chosen according to the scenarios SI 
and S2 (see Tab. |2]). For SI the grey region is excluded and for S2 the dark grey region is 
excluded. Upper left plot: tan/3 varied. Upper right plot: tan/3 = 20 and \Ai,\ varied. Lower 
left plot: tan/3 = 20 and |/x| varied. Lower right plot: tan/3 = 20 and ipA^ varied. 
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Figure 10: r(t2 — ;■ biW~^). Tree-level and full one-loop corrected partial decay widths for 
the renormalization scheme RS2. The parameters are chosen according to the scenarios SI 
and S2 (see Tab. |2]). For SI the grey region is excluded and for S2 the dark grey region is 
excluded. Upper left plot: tan/3 varied. Upper right plot: tan/3 = 20 and \Ai,\ varied. Lower 
left plot: tan/3 = 20 and |/x| varied. Lower right plot: tan/3 = 20 and ipA^ varied. 
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Figure 11: r(t2 — ;■ b2W~^). Tree-level and full one-loop corrected partial decay widths for 
the renormalization scheme RS2. The parameters are chosen according to the scenarios SI 
and S2 (see Tab. |2]). For SI the grey region is excluded and for S2 the dark grey region is 
excluded. Upper left plot: tan/3 varied. Upper right plot: tan/3 = 20 and \Ai,\ varied. Lower 
left plot: tan/3 = 20 and |/x| varied. Lower right plot: tan/3 = 20 and ipA^ varied. 
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6.2 Comparison with SQCD calculation 

Often QCD corrections to SM or MSSM processes are considered as the leading higher-order 
contributions. However, it has also been observed for SM processes (e.g. in the case of WH 
and ZH production at the Tevatron and LHC [68], for iJ + 2jet production at the LHC [69], 
or for the Higgs decay to four fermions in the SM [70]) that the electroweak (EW) corrections 
can be of similar size as the QCD corrections. Therefore, in the last step of our numerical 
evaluation, we show the size of the one-loop effects based on SUSY QCD (SQCD) only. The 
size of the SQCD corrections can then be compared to the full calculation presented in the 
previous subsection. It should be kept in mind that, following Eq. ( l44l) . also the masses of 
the scalar bottom quarks depend on the order of the calculation. Consequently, we do not 
explicitly compare SQCD with the full one-loop calculation, but analyze only the size and 
the sign of the pure SQCD corrections. 

In Fig. [12] we show the tree- level values and SQCD one-loop corrected partial decay 
widths for ^2 — )■ biH~^, t2 — )■ hiW^ , t2 — > b2H~^, t2 — )■ b2W~^, respectively. The renormaliza- 
tion scheme RS2 is used, and hard gluon radiation is taken into account. The parameters are 
chosen according to SI and S2 with tan/3 varied. For SI and S2 the grey and the dark grey 
region is excluded via LEP Higgs searches, respectively. In the lower left plot of Fig. [12] the 
curves in SI end at tan/3 f» 27 due to the closing of the phase space. The size of the SQCD 
one-loop corrections reaches the highest values for large tan/3 in the case of t2 ~^ ^1,2-f^^ 
and for intermediate tan/3 in the case of ^2 ~^ ^1,2^^^- The relative size in percent of the 
tree-level values do not exceed -8% in 1(^2 ^ biH+), +18% in 1(^2 -^ hW^), -24% in 
r(t2 — ^ b2H~^) and —6% in r(t2 — ^ b2W~^). The absolute size of the SQCD corrections 
can be compared with the upper left plots of Figs. [8]-[TT| where the full one-loop corrections 
are shown. It becomes obvious, especially in S2, that restricting an evaluation to the pure 
SQCD corrections would strongly underestimate the full one-loop corrections. (Hard photon 
radiation can be as relevant as hard gluon radiation.) Consequently, the full set of one-loop 
corrections must be taken into account to yield a reliable prediction of the scalar top quark 
decay width. 

7 Conclusions 

A scalar top quark can decay into a scalar bottom quark and a charged Higgs boson or a 
W boson if the process is kinematically allowed. These decay modes can comprise a large 
part of the total stop decay width. The decay channels with a charged Higgs boson in the 
final state form a potentially important subprocess of cascade decays which are interesting 
for the search of charged Higgs bosons at the LHC. In order to arrive at a precise prediction 
of these scalar top quark partial decay widths at least a (full) one-loop calculation has to 
be performed. In such a calculation a renormalization procedure has to be applied that 
takes into account the top quark/squark as well as the bottom quark/squark sector in the 
MSSM. These two sectors are connected via the soft SUSY-breaking mass parameter Mq 
of the superpartners of the left-handed quarks, which is the same in both sectors due to the 
SU{2)l invariance. 

Within the MSSM with complex parameters (cMSSM) we defined six different renormal- 
ization schemes for the bottom quark/squark sector, while in the top quark/squark sector 
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Figure 12: Tree-level and SQCD corrected partial decay widths for the renormalization 
scheme RS2 with tan/3 varied. The parameters are chosen according to the scenarios SI 
and S2 (see Tab. [2]). For SI the grey region is excluded and for S2 the dark grey region is 
excluded. Upper left plot: r(t2 — > biH^). Upper right plot: T{t2 — )■ hiW^). Lower left plot: 
r(t2 -^ h2H+). Lower right plot: r(t2 -^ hW^). 
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we applied a commonly used on-shell renormalization scheme, which is well suited for pro- 
cesses with external top and stop quarks. In our analysis we focused on the problem that, 
for certain parameter sets, an applied renormalization scheme might fail and cause large 
counterterm contributions that enhance the loop corrections to unphysically large values. 
We have analyzed analytically the drawbacks and shortcomings of each of the six renormal- 
ization schemes. Because of the relations between the parameters that have to be respected 
also at the one-loop level we did not find any renormalization scheme that results in reason- 
ably small counterterm contributions over all the cMSSM parameter space we have analyzed 
(we did not consider a pure DR scheme which is not well suited to describe external parti- 
cles). Some renormalization schemes (for instance, the "on-shell" scheme which is defined 
analogously to the one applied in the top quark/squark sector) fail over large parts of the 
parameter space. Others fail only in relatively small parts where, for instance, a diver- 
gence due to a vanishing denominator occurs. The most robust schemes turn out to be the 
"mb, Af, DR" (RS2) scheme and the "A;, vertex. Re Yb OS" (RS6) scheme. These renormal- 
ization schemes appear to be most suitable for higher-order corrections involving scalar top 
and bottom quarks. 

We performed a detailed numerical analysis for the full one-loop result of the partial 
decay widths corresponding to the four processes ^2 — ^ bjH~^/W~^ (j = 1,2) in our "pre- 
ferred" scheme, "?7ib, Ab DR" . The higher-order corrections, besides the full set of one-loop 
diagrams, also contain soft and hard QED and QCD radiation. We evaluated the higher- 
order predictions of the four partial decay widths as a function of tan/3, fi, Ab and ipAf We 
found mainly modest corrections at the one-loop level. Larger corrections are mostly found 
in regions of the parameter space that are disfavored by experimental constraints and/or 
charge and color breaking minima. A comparison of the full one-loop calculation with a 
pure SQCD calculation showed that the latter one can result in a very poor approximation 
of the full result and cannot be used for a reliable prediction. 

A full one-loop calculation of the corresponding branching ratios requires the calculation 
of all possible partial decay widths of the scalar top quark (and consequently a renormaliza- 
tion of the full cMSSM) and will be presented elsewhere [25] . 
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Appendix: C-functions 



As explained in Sect. 14.61 in RS6 we have to deal with infrared divergent C-functions (appear- 
ing in A(pf = 0,p^,p^)) with vanishing Gram-determinants. This case is not implemented 
in LoopTools [16]. Therefore we follow Ref. [39] (and references therein) and replace the cor- 
responding C-functions by well behaving linear combinations of B-functions q For sake of 
completeness we briefly review our implementation. The class of C-functions with only one 
external momentum zero, can be completely reduced to S-functions. Having three different 
masses we can use partial fraction decomposition: 



Co(0,p,p,mi,m2,m3) 



Bq{p, mi, ma) - Bq{p, m2, nia) 



1 1 
vn\ — VTi^ 



With only two different masses applying partial differentiation (I'Hospital) yields 

95o(p,mi,m3) 



Co(0,p,p,mi,mi,m3) 



d (ml] 



(122) 



(123) 



We also used symmetry relations and decompositions which can be found in [39] and the 
following short hand notation: 



n 4i 4i 4 0/221221 2 2\ 

Aj = P +m^ +m- - 2[p m- +p m- + m^m-) 
We included the following replacements of Cj functions with pi = 0: 

r, ,a N 5o(p,mi,m3) -5o(p,m2,m3) 
Co(0,p,p,mi,m2,m3) -^ ^ 2 ' 

mi — 7712 

Co{0,p,p, nil, 7711, 7713) -^ 7;— {p^ + ml - 77il) {2 - Bo{p, 7711,7713)) 



Ci(0,p,p,mi,m2,m3) -^ 



+ {p^ -ml- 77il)Bo{0, mi, mi) + 277ilBo{0, 7713, 7713) 

277iI{Bq{p, 7712, ms) - Bo{p, 7711,771^)) 



(124) 



(125) 



(126) 



C2{0,p,p, mi, 7712,7713) -^ 



3(m^ — mg)^ 
+ {mj - ml)Bo{p, m2, m^) -m\ + m\ 
+ (3m^ — 2m2 — 77i\ + p^)Bi{p, 7712, 7713) 

+ {ml~ml-p')Bi{p,mi,m3) 

Bi{p,mi,m3) - Bi{p,m2,m3) 

2 2 ' 

7711 ~ ^ 



(127) 
(128) 



2p 



€2(0, p,p, mi, mi, ms) -^ — -<^ 5o(0,mi,mi) -So(p,mi,m3 



p + 7111 ~ ''^3 
„2 _2 _2 



{p^ - ml + ml){2 - Bo{p,mi,m3)) 



+ {p -m^- mJi?o(0, mi, mi) + 2m^Bo{0, m^, m^) 



(129) 



^FormCalc |46j sorts the loop integrals with help of the masses. Consequently, any momentum can become 
zero, not only pi. Furthermore LoopTools uses a different convention than [39]: Ci = Cn — C12, C2 ~ C12. 
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In the case of p2 = 0, we used the following replacements: 



Co{p,0,p, mi, 1712,7713) -^ 

Co(p,0,p, nil, 7n2, 1712) -)■ 



Bo{p, mi, ms) - Bo{p, rrii, m^] 



1 



D 



12 L 



^2 — fTt^ 
I 1 2 , 2^ 

[j) — rrir.. 



Ci{p,0,p, mi, 1712, 1713) -> 



3(m 



m 



2N2 



C2ip,0,p,mi,m2,m3) -> 



2 + mi)(2 - Bo(p,mi,m2)) 
{p^ -m\- m2)i?o(0, m2, m2) + 2m^i?o(0, mi, mi) 
2ml{Bo{p, mi,m2) - Bq{p, mi, m^)) 

+ {ml - m2)5o(0, m2, mg) - m^ + m^ 
— {'im\ — 2m\ — ml — p^)Bi{p, mi, 1172) 

+ {ml-ml-p^)Bi{p,mi,m3) , 

2ml{Bo{p, mi, ma) - Bo{p, mi, m2)) 



3(m 



m 



2^2 



-iml-ml)BoiO,m2,m3) + ml-ml 
— (3m2 — 2ml ~^ '^i ^ P^)Bi{p, mi, m^] 
+ (m2 -m\- p^)Bi{p, mi, m2) 

Finally, for ps = (pi + P2) = we employed: 

ry I a N Bo{p,mi,m2)- Bo{p,m2,m3) 

Co(p,p,0,mi,m2,m3) -^ ^ 2 ' 

?7i| — m^ 

Co(p,p,0,mi,m2,mi) -^ — — {p^ -m\ + ml){2 - Bo(p,mi,m2)) 

I Jin 



Ci{p,p,0,mi,m2,m3) 



-)■ 



+ (p^ — m2 — mi)i?o(0, mi, mi) + 2mlBQ{0, 
Bo{p,m2,m3) + Bi{p,mi,m2) + Bi{p,m2,m3] 



m2,m2j 



2 2 
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D 



12 



{p^ -ml + ml){2- Bo{p,mi,m2)) 



C2{p,p,0,mi,m2,m3) -^ 



+ (p2 _ ^2 _ ^)5^(o, ^,, rni) + 2mlBo{0, m2, m2) 
2ml{Bo{p, m2, ms) - Bq{p, mi,m2)) 



3(mf — mp''2 



(2m?-m^-m^ + p2)5o(p,m2,m3) + 
— (3??ii — 2m3 — m2 + P^)Bi{p, m.2, m.3) 

+ (mj — ^^1 — P^)Bi{p, mi, m2 



2 2 

TTlg — ?Tii 



(130) 



(131) 



(132) 



(133) 

(134) 

(135) 
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